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We study a quantum dynamical system of N, O(N) symmetric, nonlinear oscillators as a 
toy model to investigate the systematics of a 1/N expansion. The closed time path (CTP) 
formalism melded with an expansion in 1/N is used to derive time evolution equations valid 
to order 1/N (next-to- leading order). The effective potential is also obtained to this order and 
its properties are elucidated. In order to compare theoretical predictions against numerical 
solutions of the time-dependent Schrodinger equation, we consider two initial conditions 
consistent with O(N) symmetry, one of them a quantum roll, the other a wave packet initially 
to one side of the potential minimum, whose center has all coordinates equal. For the case 
of the quantum roll we map out the domain of validity of the large- N expansion. We discuss 
unitarity violation in the 1 /N expansion; a well-known problem faced by moment truncation 
techniques. The 1/N results, both static and dynamic, are also compared to those given by 
the Hartree variational ansatz at given values of N. We conclude that late-time behavior, 
where nonlinear effects are significant, is not well-described by either approximation. 

PACS numbers: 11.15.Pg,11.30.Qc, 25.75.-q, 3.65.-w; Preprint number: LAUR 00-1251 



I. INTRODUCTION 

Initial value problems in quantum field theory are of 
great interest in areas such as heavy ion collisions, dy- 
namics of phase transitions, and early Universe physics. 
However, the solution of the corresponding functional 
Schrodinger equation is essentially impossible and one 
is forced to resort to approximate methods such as mean 
field approaches of the Hartree type or the large N ex- 
pansion. The application of variational techniques such 
as Hartree is limited in scope since the errors are uncon- 
trolled. While 1/N methods promise better error control 
since they are based on a systematic expansion, at next- 
to-leading order these methods can become extremely 
complicated and expensive to implement. The motiva- 
tion for our work in this paper is to implement the 1/N 
expansion at the first nontrivial order in a quantum me- 
chanical example. Not only does this simplify the analy- 
sis but it also opens the possibility of comparing the ap- 
proximate results with numerical simulations of the time 
dependent Schrodinger equation, a luxury not available 
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in the field theoretic case. However, it should be kept in 
mind that quantum mechanics and quantum field theory 
are very different. For example, in the quantum mechan- 
ics applications discussed below, the 0{1/N) corrections 
do not correspond to inter-particle collisions (as they do 
in field theory) since we are restricting ourselves to one- 
particle quantum mechanics. Nevertheless, as discussed 
in more detail below, quantum mechanical examples pro- 
vide excellent test-beds for key issues such as positivity 
violation and late-time accuracy of the approximations. 

The O(N) model has been extensively employed in 
time independent applications in statistical physics and 
quantum field theory jl]]^] and several recent applications 
have studied time-dependent phenomena. The dynamics 
of the chiral phase transition following the expansion of a 
quark-gluon plasma produced during a relativistic heavy 
ion collision has been modeled by an 0(4) a- model at 
leading order in 1/N [[^J. The nonequilibrium dynam- 
ics of an 0(A r )-symmetric Atfi 4 theory, again treated at 
leading order, has been investigated in detail Even 
at leading order, the 1/iV expansion captures the phase 
transition, but does not contain enough of the dynam- 
ics to allow for rethermalization, since direct scattering 
first occurs at next order. The O(N) model has been 
used in inflationary models of the early Universe H with 
the scalar field often starting at the top of a hill in the 
potential and "rolling" down, giving rise to a quantum 
roll problem. It has also been applied to study primor- 
dial perturbations arising from defect models of structure 
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formation j6|. 

The general method for obtaining the dynamical 1/TV 
approximation via path integral techniques in quantum 
field theory was discussed earlier in Ref. and applied 
later [^]|| to a quantum mechanical system of TV + 1 cou- 
pled oscillators (a one-dimensional version of scalar elec- 
trodynamics). Two different sets of approximate actions 
were considered, which differed by terms of order 1/TV 2 , 
both of them being energy conserving. The first method 
in Ref. || was a perturbative expansion of the gener- 
ating functional in powers of I/TV. The second method 
was to first Legendre transform the action to order 1/TV, 
and then find the equations of motion. When these two 
methods diverged from each other, they also diverged 
from an exact solution for the case N = 1. However, 
due to computational restrictions it was not possible to 
study numerically the accuracy of the approximation as 
a function of TV. Remedying that deficiency is the main 
motivation of the present study, since for the quantum 
roll problem numerical solutions can be obtained for ar- 
bitrary TV. 

One of the subtle issues in expansions involving mo- 
ment based truncation schemes such as 1/TV, which is 
present both in quantum field theory and in quantum 
mechanics, relates to the imposition of constraints aris- 
ing from the positivity of the underlying probability den- 
sity function or functional. The importance of these con- 
straints is well known in areas such as turbulence and 
beam dynamics JTo| ] . In this paper, we show that pos- 
sible violations of these constraints must be tamed in 
1/N expansions if the approximation is at all expected 
to succeed at moderate values of N. This is possible by 
using certain resummation schemes which we will discuss 
elsewhere. 

In this paper we show that the naive next-to-leading 
order 1/N expansion violates unitarity (or more gener- 
ally, positivity) leading to an instability at least for N 
less than some value TVy. We have numerical evidence 
for a sharp threshold at N ~ Nt beyond which we have 
not been able to detect an instability. This behavior ap- 
pears to be related to the nature of the effective potential 
at next-to-leading order: At this order, the effective po- 
tential has the property of not being defined everywhere 
for values of iV < N c , where N c depends on the values of 
the parameters specifying the potential, and Nt ~ TV C . 
However, for N > N c , the effective potential is defined 
everywhere. 

A comparison of the 1/N expansion and Hartree is 
of interest since both agree at infinite N. At finite TV, 
the next-to-leading order large N and Hartree approxi- 
mations differ and provide alternative routes to improv- 
ing the leading order result which, for the quantum roll 
problem, consists of harmonic oscillations in (r 2 ) where 
r is the radial degree of freedom. At finite values of N, 
the inclusion of nonlinearities leads to amplitude modu- 
lation effects on top of the harmonic motion. The ability 
to capture this modulation is a good test for the next-to- 
lcading order large N and Hartree approximations. Our 



numerical results provide evidence that neither of these 
methods are satisfactory at late times (relative to the os- 
cillation time), though they work reasonably well at short 
to intermediate times. 

Our results suggest that it is important to find ways 
to improve the naive 1/N expansion at next-to- leading 
order. Work using resummation schemes is in progress 
and short discussions of relevant issues are included in 
this paper. 

The paper is organized as follows. In Section || we 
present the O(N) mo del a s it pertains to quantum me- 
chanics and in Section III we derive equations of motion 
for the large- TV approximation to order 1/N. We de- 
rive the corresponding equations for the time dependent 
Hartree approximation (TDHA) in Section In sec- 
tion ^ we show how the same TDHA equations can be 
obtained from an equal-time Green's function approach 
which is computationally more attractive. The energies 
for the various approximations are calculated in Section 



VI. Section |VHj describes the two initial conditions which 
preserve the O(N) symmetry, namely a quantum roll, 
and the time evolution of an offset Ga ussian centered 
at an O(N) symmetric point. In Section VIII we deter- 



mine the effective potential to both order 1 /N and for the 
Hartree approximation. Numerical results and compar- 
isons with the approximations are discussed in Section 



IX 



and our conclusions are discussed in Section [x| 



II. THE O(N) MODEL 

The Lagrangian for the O(N) model in quantum me- 
chanics is given by: 
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i=l 



where V(x) is a potential of the form 



(2.1) 



V(r) = -?- (r 2 -r 2 ) 2 , 
W 8TV y °> ' 



N 
= 1 



X: 



(2.2) 



The time-dependent Schrodinger equation for this prob- 
lem is given by: 



dip(x, t) 
di 



f i N a 2 1 

p|X^ +V(r)^{x,t). (2.3) 



For arbitrary initial conditions, given present computa- 
tional constraints, these equations can be numerically in- 
tegrated only for small N < 4. The initial conditions for 
the quantum roll problem allow a numerical solution for 
all N, and in this case we can attempt to study fully 
the behavior of the large- N expansion. (For the shifted 
Gaussian initial conditions, however, this is not possi- 
ble, and we used numerical solutions obtained for TV = 1 
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and 2 to benchmark the large- N approximations and the 
TDHA solutions at short times.) 

The symmetry of the quantum roll problem is such that 
only the radial part of the wave function is of interest. 
Assuming a solution of the form 



ip(r,t) = r( 1 - N V 2 (f>(r,t), 



(2.4) 



the time dependent Schrodinger equation for <fi(r, t) re- 
duces to plfl : 



dt 



with an effective one dimensional potential U (r) given by 
(N-l)(N-3) , g 



U{r) 



8r 2 



It is further useful to make the rescaling: 

r 2 = Ny 2 , r 2 = Ny 2 . 
The potential (|2.6|) then becomes: 



(2.7) 



f/(y) (JV-lKjV-3) o 2 2,2 



N 



8N 2 y 2 



(2.8) 



corresponding to the new Schrodinger equation, 



dt 



i d 2 

' 2N 2 dy 2 



+ u(y,N)\ <t>(y,t) (2.9) 



where t — N t. 

The method of choice to investigate the long-time be- 
havior of the exact solution is the split-operator method, 
which has been presented in detail in Ref. JT2J . The wave 
function is expanded as a Fourier series in the radial com- 
ponent, and the solution is obtained as the repeated ap- 
plication of a time-evolution operator in symmetrically 
split form. As a result, the use of a Fast-Fourier Trans- 
form algorithm is required. For the purpose of the present 
implementation, 256 radial grid points, a value of 20 for 
the radial grid boundary, and a time step size of 0.01, 
provides a conservation of the wave function unitarity to 
better than 9 significant figures. The accuracy of the 
method has been established by comparing results with 
a second method, where we first solve for the eigenvalues 
and eigenfunctions, and then use the expansion 



E 



C n e 



-%E„t 



(2.10) 



where C n was determined from the initial conditions. 
This method was restricted to moderate values of N. 
Results from the two methods agreed in the cases where 
they were used together. 



III. THE LARGE- TV APPROXIMATION 

The large N approximation has been worked out for 
the O(N) mode l in 1 + 3 dimensions in Ref. p]. The 
Lagrangian (2.1) with the potential function ( |2/2 ) is ob- 
tained from that paper by specializing to + 1 dimen- 
sions, and replacing <f> a {t) ~^ Xi{t), v — > Vq, and A —* g. 

To implement the large N expansion, it is useful Q to 
rewrite the Lagrangian in terms of the composite field \ 
by adding a constraint term to (p-l|) , given by: 



N 
2f) 



2N 



(r 2 r 2 ) 



which yields an equivalent Lagrangian, 



L'(x,x,x) = 



2\ r N 2 



(3.1) 



(3.2) 



The generating function Z[j, J] is given by the path in- 
tegral over the classical fields Xi(t): 



Z[j, J] 



JW[j,J] 



dx Yl dx * ex p{* s i x > x; 3, J] } j 

S[x,x;j,J]=J^ dt |i' + 3iXi + Jx} • 

The effective action, to order 1/N, is obtained by inte- 
grating the path integral for the generating functional 
for the Lagrangian (3.2), over the Xi variables, and ap- 
proximating the integral over x by the method of steep- 
est descent (keeping terms up to order 1/N). A Legen- 
dre transform of the resulting generating functional then 
yields the effective action, which we find to be: 

r[g,x] = 

/ [#(*)- X(*)9?(*)] +~E ln [G« 1 (V)] 

Jc I „■ 



r N i 

-fx(t) + ^x 2 (t) + ^HD-%t)} 



(3.3) 



where the integral is over the close time path C, dis- 
cussed in Ref. [| and q(t) = (xi(t)). Here G^(t, t') and 
D^ 1 (t, t') are the lowest order in 1 /N inverse propagators 
for Xi and x, given by 

G^^t') = |^ + X (t) j S c (t,t')6 l3 = G-\t,t')6 l3 , 

D- 1 (t,t r ) = --8 c (t,t')-U(t,t') : 
9 



where 



U(t,t') 



-Y,G l3 {t,t')G Ti (t',t) 
Qi(t) G tJ (t,t') qj (t'). 



(3.4) 
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Here Sc(t, t') is the closed time path delta function. 

The equations of motion for the classical fields qi(t), 
to order 1/N, are 



dt 2 



+ X(t)>q t {t) 



+i V f dt'G ij (t,t')D(t,t')q j (t')=0, 



with the gap equation for x(t) given by 



x(t) 



9 

' 2N 



9 

2N 



E 



(2) 

The next-to-leading order Xi propagator Q\- (t,t ) and 
self energy Ey(t, t') to order 1/N turn out to be 



(3.5) 



(3.6) 



d 2) (M') = G«(t,f) 



(3.7) 



-J2 I d< i / d*2G <fc (t,ti)E w (ti,«2)G w (t 2 ,t / ), 

Ejuftt') = iG H {t,l/)D{tJ)-q k {t)D{t,l/)qi{l/). 

These equations agree with (2.18-2.22) of Ref. @. We 
mention here that the actual equation for Q w hich follows 
from the effective action differs from Eq. ( |3.7|) in that the 
final G in the integral equation is replaced by the full Q. 
This leads to a partial resummation of the 1 /N correc- 
tions which which guarantees positivity of (x 2 (t)) (this 
restricted result does not imply that the full positivity 
problem for the density matrix has been solved). How- 
ever, it does not improve the long-time accuracy of the 
results HI- 

In order to solve for D(t,t'), we first write 

— D(t,t') = -S c (t,t') + — AD(t,t') . (3.8) 
9 9 

Then AD(t,t') satisfies the integral equation, 

-AD(t,t') = ^n(M')- / dt"U(t,t")AD(t",t'), 
9 N J c 

(3.9) 

in agreement with (2.13-2.16) of Ref. @. 

We are now in a position to solve these coupled equa- 
tions for the motion of qi(t) and x(£) for given initial con- 



ditions. For the initial conditions discussed in Sec. [VTj, 
we find: 

G ij (t,t')/i = 6 ij f(t)f*(t'), (3.10) 
where f(t) and f*(t) satisfy the homogeneous equation, 



d 2 

dT 2+ ^ 
with initial conditions: 



/(*) 
/*(*) 



0, 



(3.11) 



/(0) 



G , 



/(0) = 1/(2VG) 



(3.12) 



However, AD(t,t') cannot be factored into products of 
functions like G%j (t, t' 



W e solve ( |j.5| ) and (3.6) simultaneously with (|3~7|) and 
( |3.9|) , using the Chebyshev expansion technique |l3|] of 
Appendices A and B of Ref. B . 



IV. THE TIME DEPENDENT HARTREE 
APPROXIMATION 

It is useful to compare our results for the large- N ap- 
proximation to the time dependent Hartree approxima- 
tion (TDHA) suitably formulated for the O(N) problem. 
The static Hartree approximation is based on the idea 
of varying the parameters of a Gaussian wave function 
so as to minimize the energy (the generalization to the 
time-dependent case is given below). For the O(N) prob- 
lem this amounts to placing an iV-dimensional Gaussian 
some (radial) distance away from the origin and then car- 
rying out the minimization procedure. In contrast, the 
leading-order large N wave function is a Gaussian which 
is locked at the origin. At infinite N the TDHA becomes 
exact and equivalent to the leading order large- N approx- 
imation, a well-known result. At finite N, the TDHA and 
the next-to-leading order large N approximation can be 
thought of as two competing schemes to improve on the 
leading-order result. 

There are several ways of implementing the Hartree 
approximation: The most common is by using the time- 
dependent variational principle of Dirac p4[~^6[. This 
has the advantage of giving a classical Hamiltonian de- 
scription for the dynamics of the variational parameters, 
which can be hidden in other formulations. 

The idea behind this approach is that the variation of 



dt(ip(t)\i§- t -H\m) (4-1) 



is stationary for the exact solution of the Schrodinger 
equation. We consider Gaussian trial wave functions of 
the form: 



ip{x, t) = Af cxp 



ipi(t) Zi(t) 



(4.2) 



-Zi(t) 



-Hy(t) z 3 {t) 



where Af is the normalization constant, and we have set 
Zi(t) = Xi - q l (t). Here qi(t), pi(t), G lj (t) and Hy(t) are 
time-dependent variational parameters, to be determined 
by minimizing the Dirac action. We note that Hij(t), 
which is used only in this section, is conjugate to Gij(t) 
and is not to be c onfused with the self energy U(t,t') 
defined in Eq. fljuj) . 

The n-point functions can be calculated from the gen- 
crating functional using the formula, 
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(ZiZj ■ • • z n 



d n Z[j] 



djidjj ■ ■ ■ dj„ 



(4.3) 



3=0 



where 



Z[j]=M 2 



J]_dx s 



x exp 



-z^G^z^t) + j iZi (t) 



exp 



jiGijjj 



(4.4) 



The expectation value of the time derivative is given by 



d 



and the expectation value of the kinetic energy is 
1 d 2 \ ViVi 1 



2 dx 2 



Gl 



2 TlijGjkTlki ■ 



(4.5) 



(4.6) 



For the expectation value of V we first expand the po- 
tential in a Taylor series about Zi = 0, 



V{q, z) = V(q) + Vi(q) Zl + -Vij(q) z iZj + 



where 



2N 
9 

2N 



r. 



[Sij (q s q s 



o) > 

'2\ 



(4.7) 



r, 



o) 



V i3 k{q) = (<>i3<ik + $ik<lj + Sjkli) , 
Vijki{q) — jq {dijSki + SuSjk + SikSji) . 

Thus 

V{q, z) = jyy[(qjqj - rlf + (z lZl ) 2 + 4{ziZi){zjqj) 
+ i(ziqi) 2 + 2{z i z i ){q j q j - Tq) + ^{mililjlj - r%)] 
Taking the expectation value, we obtain: 



(V) = ^[{QiQi ~ r oV + 2G u (q jqj - r 2 ) 
+ AGijqiqj + GaGjj + 2GijGj%\ 



(4.8) 



The Hartree equations of motion are Hamilton's equa- 
tions for the variational parameters: 



Qi=Pi, 
Pi = —Vi 



(4.9) 



1, 



-VijkGjk , 



Gij — 2 (GikHkj + GjfJIki) > 



Vijki Gu 



Solutions of this set of equations determine the time- 
dependent Hartree approximation to the true solution of 
the Schrodinger equation. 



V. METHOD OF EQUAL-TIME GREEN'S 
FUNCTIONS 



Solutions of the TDHA equation (|4.9| ) require com- 
puting the matrix inverse of Gy . This can be difficult 
to carry out in practice for large N. Fortunately, the 
method of equal-time Green's functions provides a way 
to avoid this technical difficulty (l7]Jlj| • We begin by con- 
sidering the time-evolution of the one point functions: 



Qi = V% , 

Pi = —Vi — — Vijk Gjk 



(5.1) 



2N 



{qiiqkqk + Gkk - to) + qk(G ik + G ki )} 



where Vi and Vijk are given by (4.7), as well as the evo- 
lution of the two-point functions: 



G ij(t) — (zi Zj) , Kij — 
F ij = \{ \ Z i*3 + *3 Z i\) ■ 



(5.2) 



Here we have again set Z{(t) = x, — qi(t). All of the 
expectation values are tak en w ith respect to the Gaussian 
trial wave function, Eq. (4.2). To obtain the equations 
of motion for the two-point functions, we use the exact 
equation of motion and the factorization resulting from 
the Gaussian approximation. This yields 



Gij — Fij 



F 



Ki 



dV \ 
dxj I 



(5.3) 
(5.4) 



Kij = 



[ dV 

dx,. 



dV 

dxj 



(5.5) 



Here we have used the Lagrange equations of motion 
dV 



dx, 



= 0. 



(5.6) 



where 



9V T7 1 T T 

— Vij Zj + - Vtjki Zj z k zi 

+ terms with even powers of z^ , 



(5.7) 



and the fact that for our Gaussian wave packet, (z-i) = 0. 
The canonical commutation relations give: 

(ZjZj ZjZi) — (xiXj XjXi^j — iSij , (5.8) 

and we have 

{zi Zj z k zi) = GijGki + GuGjk + GikGji , (5.9) 
(zi Zj Zk zi) = GijFki + FuGjk + GikFji 

+ i {GijSki + GjkSu + dkSji) 12 , (5.10) 
{Zi Zj Zk zi) — FjiGki + FuGjk + FkiGji 

- i (GkiSji + GjkSu + GjiSki) /2 . (5-11) 
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Finally, from Eqs. (5.4) and (|5.5|) we get: 



F; 



Ki 



VjkGik 



(5.12) 

— Vjklm {GikGim + Gi m G k l + GuGkm) /6 , 

Kij = -V ik F kj - V jk F ki (5.13) 

— Viklm {GklFmj + GlmFfoj + GkmFlj) /6 
Vjklm (FkiGlrn + F m iGkl + FuGkm) /6 . 

For Gaussian initial conditions, the equal-time Green's 
function method is assured to give the same result as the 
Hartree method, if i*#(0) and ifjj(O) satisfy the require- 
ments: 



F«(0) = ; 



^(0) = Gr>(0)/4 



Choosing .KV,- independently of G^, corresponds to a 
mixed initial density matrix, rather than a pure state. 
If we further choose Gjj(O) to be diagonal and equal to 
the same number Go, 

Gij (0) = Sij Gq , 



then Kij(0) is given by: 



K l3 (0) = 6^/(4 G ). 

For the initial conditions pertinent to the quantum roll, 
q l {t) = for all t. Then G lj (t), F %j (t), and K tj {t) are all 
proportional to the unit matrix, and have no off-diagonal 
terms. For the offset initial condition, we choose large- 
TV symmetric initial conditions so that qi(0) = qo, and 
Pi(0) = 0, and Gy(0) = Go<5y. In that case, all the q's 
and p's are identical, 

q t (t)=q(t), Pi(t)=p(t), 

and the matrices Gij(t), Fij(t) and Kij(t) become off- 
diagonal in a simple way so that all the diagonal el- 
ements are equal and all the off-diagonal elements are 
equal. That is, we can write 

G ij (t) = G(t)8i j + G(t)Q.-6 ij ), 
F ij (t)=F(t)5 ij + F(t)(l-5 ij ), 
K ij {t)=K{t)5 ij + K{t){l-8 ij ). 



For this case, Eqs. (fT|), (0), (HI), and |yj) simplify 
to the following set of coupled equations: 



Q = P, 



(5.14) 



9_ 

2JV 

2F, 



q<Nq 2 - rg + (JV + 2)G + 2(N - 1)G 



G = 2F . 



G 

F = K ^[g[{N + 2){q 2 + G) -r 2 \ 



+ 2{N -l)G(q 2 + G) 



F = K - ^- |g [ (37V - 2) g 2 + (JV + 4) G 



2(JV-2)G-rg] +2Gg 5 



X = --^F[(JV + 2)(<? 2 + G)-r 2 l ] + 



+ F2(7V-l)( g 2 + G) 



if = --^<j F [ (3JV - 2)g 2 + (JV + 2) G 



+ 2(JV-2)G-r 2 ] + F2(g 2 + G) 



If we let rg = JVt/o) and then take the limit JV — > oo, 
we recover the leading order in the large- JV result, as 
discussed in Ref . |jl7) . 

The equal time Green's function method is easier to 
implement numer ical ly than the Hamiltonian system de- 
scribed by Eq. (4.E), since no matrix inversion is in- 
volved. However, if one wants to find the wave function 
or the energy, instead of just obtaining the Green's func- 
tions, matrix inversion is once again required. 



VI. ENERGY 

It is important to note that even though the Hartree 
and large JV approximations are truncations of the true 
dynamics, they are nevertheless energy conserving. In 
the large- JV approximation, to order 1/JV, the expecta- 
tion value of the Hamiltonian is given by 

E = ^£[te a (*)> + <*(*) *?(*)>] (6-1) 



5 (*(*)> 



JV 
2^ 



(r(t)) 



We write these expectation values in terms of the CTP 
Green's functions. By definition, the disconnected two- 
point Green's functions are introduced as 

D dia (t,t')=i(T e [x(t)x(t')]) 

= i X (t)x(t')+V(t,t') 7 (6.2) 

Gij t dis(t,t') = i{Tc[xi(t)xj(t')]) 

= iq i {t)q j (t')+g ij {t,t'), (6.3) 

where V and Q a b denote the connected two-point Green's 
functions 



V(t,t') = 



S 2 W[J,j] 
8J{t) 8J(t') 

8 2 W[J,j] 



J, 3=0 



J, 3=0 



(6.4) 



To obtain the energy from ( J6.l| ) , we require the expecta- 
tion values 
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(x(t)) = x(t) , 

(x 2 (t)) =q 2 (t) + Gu(t,t)/i 



(x 2 (t))=q 2 (t) + 



dtdt 1 



(X(t)x 2 (t)) =x(t) [q 2 b {t) + gu{t,t)/i] -Kii(t,t,t), 

where Kij(ti,t2,ta) is the 3-point Green's function de- 
fined as 



Kij(ti,t2,t3 



dtG ik {h,t) G kj {t 2 ,t) D(t,t 3 ) 



The energy for the next-to-leading order large-N approx- 
imation is then given by: 



E 



N 



° T x(t)- 7r Ax 2 (t) + AD(t,t)/ l } 



2fj 



d 2 Gu(t,t')/i 



dt dV 



(6.5) 



Using the equations of motion, we can show directly that 
( |6.5D is conserved. 

It is easy to evaluate the energy at t — fo r the quan- 
tum roll problem using the initial conditions ( |3.12 ) . The 
result for the leading and next-to-leading order large-N 
approximation is: 



where, 



E _ 



8G + 8 9Vo 

1 r 2 



1 

iv £l 



- A 9Gyl 



(6.6) 



Our initial wave function was chosen to be Gaussian, so 
that the parameters of the Hartree approximation agree 
exactly with the energy and parameters of the exact wave 
function at t — 0. However, the leading order in the large- 
N approximation for the same value of G will disagree 
with the exact energy by ei . This discrepency dissapears 
when we include the 1/N corrections. 

Since the Hartree approximation leads to a canonical 
Hamiltonian dynamical system, the corresponding en- 
ergy in that approximation is also a constant of the mo- 
tion. It is given by: 



-G- 



+ 2HijGjk^lki 



(6.7) 



+^jVyfcz {GijGki + GuGjk + GikGji) 



We used this expression to check the accuracy of our 
numerical solutions. As with the next-to-leading order 
1/N exp ression, for the quantum roll initial condition, 
Eq. (|6.7[) agrees with the exact result. 



VII. INITIAL CONDITIONS 

A. Quantum Roll 

We wish to study initial conditions which are consis- 
tent with O(N) symmetry. This implies immediately 
that all the Xi(t) have to be identical, with Xi(0) = 0, 
and Gij (t) must be diagonal. The quantum roll prob- 
lem is defined by a Gaussian initial wave function that is 
centered on the origin: 



i>o{r) 



1 



^TrG)^/ 4 



exp 



r 
~4G 



(7.1) 



In this section, G = G(0). 

One of the difficulties in studying the systematics of 
the 1/N expansion is the fact that, at next to leading or- 
der, every different value of N (with all other parameters 
held constant) defines a different initial value problem. 
In this sense one cannot naively compare individual so- 
lutions, exact or approximate, at different values of N. 
In effect one has to tune the parameters of the problem 
at each TV in order to maintain certain invariance proper- 
ties which allow different N evolutions to be compared to 
each other. This parameter tuning process is described 
below. 

Since the infinite N limit has very precise properties, 
several technical issues arise when one wants to approach 
this limit starting at N = 1 in a uniform manner. To 
study the large TV limit it is conv enie nt to make a rescal- 
ing to the y variables, given in (|2.7| ). At very large N, 
the potential energy u(y,N) is (|2.S|): 



u{y,N) 



(JV-l)(iV-3) , g 



2\2 



8/V 2 J/ 2 



+ £ (v 2 - vt) 



2\2 



(7.2) 



In this limit, u(N,y) has a minimum which is indepen- 
dent of N, and the large N limit consists of harmonic 
oscillations about this minimum [the reason for this is 
that the large N limit also corresponds to an effectively 
large mass limit in the Schrodinger equation ( |2.9| )]. 

One way to uniformly study the motion of a wave 
packet as a function of N is to choose initial conditions 
so that there is a uniform overlap of the initial wave func- 
tion with the set of eigenfunctions of the Hamiltonian in 
the N — > 00 limit. We can obtain this constant overlap if 
we allow the coupling constant g to be a slowly varying 
function of N. This can be done in several ways that dif- 
fer by terms of order 1/N 2 . The method presented below 
leads to uniform results even at N = 1 as we change the 
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parameters with A. Our method is to keep the distance 
between the centers of the initial wave function and the 
position of the minimum of the potential a constant as 
A is varied. 



Using (7.1), we define f by: 

T((tf + l)/2) 



'2G 



and Gqo by the variance: 



T(N/2) 



(7.3) 



= r~ - r 



= G^ A-2 



r((A+l)/2) 

r(iv/2) 



Solving the above equations for G, we have 

Goo 



G(A) 



2N - 4 [T((N + l)/2)/r(JV/2)] 
Substitution of this expression into ( [7.3[ ) yields 

f(JV) 



G 



AT[r(JV/2)/r((iV+ l)/2)] 2 - 2 
In the limit when A goes to infinity, we have 
G(A) 



(7.4) 



(7.5) 



G 



f (AO = - 1)G 



which defines foo, and agrees with the asymptotic form 
of the rescaled version of the initial wave function ( |7.1| ) : 



Mr) 



1 



1 

(2ttG) n I a 



exp 
exp 



AG 



A- 1 



lnr 



(r - Too) 2 
2Goo 



0(1/VA) 



In order to ensure that the initial wave function has 
a finite overlap with the energy eigenfunctions of the 
Schrodinger equation at large N, we will keep the value 
of Goo (and not G) fixed in our simulations. 

Another quantity that should be kept constant is the 
basic oscillation frequency. In order to do this, we first 
find the Gaussian oscillations about the mini mum of the 
one dimensional potential, defined by Eq. (f?~6|). Wc 
expand U(r) as 



U(r) = U{f) + im 2 {r-ff + ■■■ , 

where f is given by the solution of the equation, 
(A-l)(A-3) g 



4f 4 

and m 2 by 

3(A-l)(A-3) 



2N 



(7.6) 



(7.7) 



9 /„ 

]v (3r 



4f4 

2 o „2\ 



+ ^7(3r 2 -r 2 ), 



2N 



(7. 



24) 



The frequency of oscillation is determined by m, and this 
is the quantity to be kept fixed as is changed. 

The last technical issue is to keep the distance between 
the center of the initial wave function, f, and the mini- 
mum of the potential, f, a constant as we vary N. That 
is, we keep 



5r 



constant for all N. 

With this strategy of keeping Goo , fh and Sr fixed, we 
can now determine how the coupling constant must vary 
with N. We first define m 2 to be the second derivative 

f, 



of V(r) evaluated at r 
^2 d 2 F(r 



dr 2 

Then, ( [7.8[ ) becomes 

m 2 (N) = to 2 - 



9 

2N 



(3f 2 



3(iV-l)(JV-3) 



4f 4 



(7.9) 



(7.10) 



Solving (7.9) for ro, substituting into ( |7.7| ) and solving 
for g gives: 



(A — 1)(N — 3) 



(7.11) 



with r = f + Sr. The value of r 2 is then determined by 
_ _ 2 2A (A-l)(A-3) 



(LZ) 



4f 4 



(7.12) 



Thu s, for fixed v alues of Goo, m and Sr, Eqs. ( [7.4| ), 
Q , ( PH) , and ( [ml ) determine values for G(A), f (A), 
5(A), and r (A) all values of A. 
In the limit, A — * oo, wc find that 



m 2 (oo) = m 2 



4G 2 



(7.13) 
(7.14) 



To summarize, in order to establish appropriate initial 
conditions for the quantum roll problem, we have kept 
the variance Goo constant instead of G, and have allowed 
the parameters describing the potential function g and ro 
to change with A in order to compare solutions that have 
close to the same oscillation frequencies. In our numerical 
runs, we chose the values 



Goo = 1; m 2 = 2; Sr = 2; g(oo) = 1 



(7.15) 



Fig. [I] displays the variation of the potential parameters 
with A. 
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FIG. 1. Potential parameters as a function of N. 



B. Shifted Gaussian Initial Conditions 

The second O(N) invariant initial condition we inves- 
tigated had a wave function localized in a wave packet 
near the center of the valley of the classical potential at 
r = ro. For N = 1 this would be the standard double- 
well tunnelling problem; for higher values of N, tunnel- 
ing is avoided by going around the barrier. Therefore 
this initial condition is qualitatively different from the 
roll problem and provides a different arena for testing 
approximations. However, since this initial condition vi- 
olates the O(N) symmetry of the potential, numerical 
solution is at present possible only for very small N. 

We take the initial wave function to be a shifted Gaus- 
sian of the form: 



ip (x) 



(2nG) N / i 6XP { ^ 



{Xj - r /VN) 
AG 



(7.16) 



The ene rgy E of this state can be determined from 
Eq. (6.7) by the substitutions, 



G, 



' 'J T ^ij ^ i 

from which we find: 

E = ^ + ^{ N ( N + 2 )G 2 -2NGrl+rtY (7.17) 

On the other hand, the height of the classical potential 
barrier is given by: 



r /VN; pi 







n, 



o, 



TP 9 4 
Eh = Tn . 

87V 



(7.18) 



For N = 1, the necessary requirement for tunneling is 
that E < E b . 



In the general case (arbitrary N), we have: 



M 



d 2 V(r) 



dr 2 



N 



N 
fJ 



M 2 . (7.19) 



If the initial state is close to the ground state of a har- 
monic potential that approximates the potential at the 
bottom of the well, then the width G of the wave function 
is, approximately, 



G 



1 



(7.20) 



which can be combined with Eq. (7.17) to give the desired 
energy of the initial state in terms of the values of N and 

a- 

We are interested in initial conditions where the energy 
per oscillator does not increase as a function of N. To 
implement this we fix M 2 = 1, which corresponds to 
G = 1/2 for the initial width. The barrier height is then 
given by 



and the total energy by 



E = 



N+l N + 2 
4 + 32 



-9 



(7.21) 



(7.22) 



We explored three cases: E — 0.5 Eh, E = Eh , and 
E = 2 Eb- For each of these cases, Eqs. ( 7.21 ) and 
( |7.22[ ) determine g for each N. In all cases we took 



Xi(0) - r /VN, ±i - 0, G tf (0) = GSij, and G« (0) = 0. 
As a consequence, all of the oscillators Xi(t) move iden- 
tically. 



VIII. EFFECTIVE POTENTIAL 

It is well-known that the static effective potential is not 
always a useful guide to the true dynamics of the system 
(See, e.g., Ref. [Q). Nevertheless, one may seek to gain 
qualitative insight into some aspects of quantum dynam- 
ics this way, though care is certainly indicated (See, e.g., 
Ref. jl9| for the Gaussian effective potential). Indeed, 
there appears to be an interesting connection with the 
properties of the effective potential at next-to-leading 
order and with the corresponding dynamical evolution 
(discussed in the next section). 

The effective potential in the large N approxima- 
tion has been previously obtained by Root [^0) to order 
1/N; however, we recalculate it here using our equations. 
When Xi and \ are independent of time, we can ignore 
the closed time path ordering and use Fourier transforms, 
passing the poles by using the Feynman contour. Then, 
from the action given in (|3.3|), we find 



9 



XT 



(8.1) 



where % satisfies the requirement 
(9 



dx 



V eS {r, X ) = 0. 



(8.2) 



In this section to make contact with Root |2C|| . we have 
/i 2 = —gr^/(2N) < 0. In o rder to examine the large AT 
limit, we again rescale (2.7) to the y variables. Then for 
the Green's functions, we find: 



There exist two real solutions of Eq. (8J3) for \ with 
y greater than some minimum value y m i n - The nex t-to - 
leading order large N effective potential, from Eq. (8.4) 
is therefore double valued for y > y m i n , and does not ex- 
ist for smaller values of y. The physical solution branch 
corresponds to the one that matches on to the leading 
order result; the other branch is an unphysical solution. 
Since it follows from a Legendre transformation, the ef- 
fective potential (at any order in 1 /N) has to be a convex 
function. The nonexistence of the effective potential at 
y < y m i n implies that no quantum state can be associated 
with the next-to-leading order large N approximation in 
this range. 



G- 1 (k) = -(k 2 - x ) 
N 



iN 



dp p. 



D-\k) = N y* G {k) + — / ^G(p)G(k-p) 

g 2 J 2n 

= -Eli- y 1 9 1 ) 

9\ 9 k 2 -x 2Vxk* -4 X J 
N (k 2 - m\){k 2 - mV) 
= ~J (k 2 - 4 X )(fc 2 - X) ' 

where mi. = b ± \/b 2 — c, with 



1 



2VX 
1 



c = ^x +g[^y x + ^Vx 

For the Feynman contour, we have 

[ ln(fc 2 — x) = Vx + constant terms . 
J 2?n 



Thus the effective potential (3.1) becomes 



- 2 [y y ) 2g + 2 



N 



(8.3) 



(8.4) 



1 

2N 



(m+ + m_ — 3 y/x) ■ 
The gap equation which determines x follows from ( |8.2| ) 



9(2 2\ _, 9(N-3) g d(m++m-) 
X =2( y - J ^ + ^N^ + 2N d~ X • (8 - 5) 



To leading order in the large N expansion, Eqs. (8.4, 8.5) 
reduce to the parametric set 



v^ix) _x 2 Vx 

N 2g 4 
2 

y 2 (x) =y 2 + -x- 
g 



2VX 



(8.6) 



Equations (^J) and ( |8.q ) agree with Root, however he 
used the leading o rde r expression for x m ( ^-6[ ), rather 
than the full x of (|J). 
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FIG. 2. Vcs/N vs y — r/yN for the leading and next-to- 
leading order large-iV approximation for different values of 
AT. 

In Fig. ||, we plot the physical branch of the effective 
potential as a function of y, for values of N from 1 to 
100, for the case g = 1 and yo — 2. For comparison, 
we also show in t his fi gure the leading order potential 
function from Eq. (S.6), which is single-valued and finite 
for all y. (In contrast to the next-to-leading order case 
we can always associate a Gaussian wave function with 
the leading order approximation.) 

In the case of the Hartree approximation, one can de- 
fine an "effective potential" as the expectation val ue o f 
the Hamiltonian using the variational wave function ( |4.2| ) 
for static configurations |l9|j2l| ]. Setting Pi(t) = and 

,qf=r 2 and Gy = #y G in 



IL, (t) = 0, and putting 



Eqs. (4.6) and (4.8), we find 



v£ ] (y,G) 

N 



1 

8G 

+ g- 



g ( 2 
~ 8 (y 

N + 2 
47V 



2\ 2 

Vo) 



(8.7) 



G G 



vlG. 



The value of G is fixed by the requirement that 



10 



dG 



0, 



which gives the gap equation for the Hartree approxima- 
tion: 



X 



9(2 2n , 9 I 2 , 1 



fj 



4VX 



where we have set G = 1/2^/x- Parametric equations for 
the Hartree effective potential are then given by: 



v [ l\x) 

N 



1 



N 



2g \N + 2 

N + 4 
+ 4(iV + 2) 

9 



N 



(N + 2) 



2 VoX 



Vx-- 



2(N + 2) 

9_i 

4(iV + 2) y/x 16 7V X 



2 VO 



(8.9) 



y 2 (x) 



N 



TV- 



1 



2 V %+ ~9 X 



2VXJ (N + 2)x 



Note tha t in the limit N — > oo, Eq. ( [3.9] ) reduces to 
Eq. (8.6), the leading order large iV result. Note also 
that the Hartree effective potential is not derived from 
a Legendre transform and hence is not subject to a con- 
vexity constraint. 
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FIG. 3. Vefi/N vs y = r/y/N for t he Ha rtree approximation 
using the parameters found in Eq. (7.15). 



In Fig. ||, we plot the Hartree effective potential from 
Eq. (3.9) as a function of y, for our chosen parameters 
goo = 1 and yo — 2, for different values of N. In contrast 
to the smooth behavior exhibited by the large N effective 
potentials, the Hartree effective potential shows a "first- 
order transition" in the placement of the minimum of the 
potential as a function of N. 

The minimum of the effective potential corresponds to 
a determination of the ground state energy. In Fig. ^, 



we show values of the minimum energies of the large N 
and Hartree effective potentials as a function of N. For 
N > 2, the Hartree minimum is generally greater than 
that for the next-to-leading order large- N approximation. 
Since the Hartree approximation is a variational ansatz it 
gives an upper bound to the minimum energy. The fact 
that the next-to-leading order large N results are lower 
than this bound is encouraging, although no guarantee 
of absolute accuracy. 
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FIG. 4. V e s/N at the minimum for the Hartree and 2nd 
order large- iV effective potentials, as a function of N, for the 
same set of parameters as in Fig. H. 
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FIG. 5. The value of y at the minimum of the effective po- 
tential for the Hartree and 2nd order large- iV approximation, 
as a function of N, for the same set of parameters as in Fig. []. 

It is interesting to ask the question how the point y m - ln , 
below which the next-to-leading order large N effective 
potential does not exist, changes as a function of N. We 
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know that at "infinite TV," (leading order), y m i n = 
but it is important to know how this limit is reached. 
For instance, is there a finite value of ./V beyond which 
y m in = 0? In Fig. ||, we plot j/ m ; n as a function of N, for 
the Hartree and next-to-leading order large TV approxi- 
mations. As already stated, the Hartree approximation 
displays a first order phase transition between the broken 
and unbroken symmetry solutions at TV — 6.2, whereas 
for the next-to-leading order large TV approximation a 
different type of behavior is found: for TV < 18.6, y in ; n 
is finite, but for TV > 18.6, it hits the origin. Thus for 
TV > 18.6, we can associate a quantum state (though not 
known explicitly) with the next-to-leading order approx- 
imation. 

The critical value of TV is fixed by the value of x at the 
gap equat ion at the inflection point. If we write the gap 
equation (jglj) as 



f(x,y 2 ,N) = f Q ( x ,y 2 ) + ^fi(x,y 2 



where 



Mx,y 2 ) 
h(x,y 2 ) 



■ (y 2 - vl) + 

3g 1 



4 VX 
g d{m + 
2 







X , 



(8.10) 



9x 



Then the critical point is determined by 

h(x,o) dh( x ,o)/dx 



fo(x,0) 9/ (x,0)/9x' 



(8.11) 



where x is given by the solution o f this system of equa- 
tions. For the parameters of Eq. (7.15), we find numer- 
ically that TV C — 18.60 in excellent agreement with the 
results shown in Fig. [|. 



IX. NUMERICAL RESULTS 

A. Quantum Roll 

We begin with a discussion of our results for the quan- 
tum roll problem. We first examine the short time be- 
havior, < t < 3, to see if the next-to-leading order 
large TV approximation gives an improvement over the 
leading order solutions. In Fig. ^, we plot the values 
of (r 2 )/N from the numerical solution, the leading and 
next-to-leading order large- N approximations, and the 
Hartree approximation, for N = 20. The next-to-leading 
order large N approximation is clearly better than the 
leading order solution, and also better than the Hartree 
results. Similar behavior is seen for other values of N 
(we also ran N = 50, 80, and 100). 
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FIG. 6. (r 2 )/N for TV = 20, for the exact, leading and 
next-to-leading order large-TV, and Hartree approximations 
for short times. 
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FIG. 7. (r )/N for the exact, leading and next-to- leading 
order large-TV, and Hartree approximations for long times. 
The top figure is for TV = 20 while the bottom figure is for 
TV = 100. The labeling conventions are the same as in the 
previous figure. 

The long time behavior of these approximations is typ- 
ically of much more interest. We examined behavior over 
times < t < 100 to see how long the approximations 
remained viable. Fig. ^ displays (r 2 )/N for the numeri- 
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cal solution, the leading and next-to-leading order large- 
TV approximations, and the Hartree approximation for 
N = 20 and 100. The next-to-leading order large-TV 
approximation for TV = 20 blows up at t ~ 84. This 
instability is connected to a violation of unitarity in the 
particular implementation of this approximation and will 
be discussed in greater detail below. In general, at these 
moderate values of TV, the approximations track the nu- 
merical solutions reasonably well though they do get out 
of phase as time progresses. As TV is increased, the phase 
errors are considerably reduced as is apparent in the re- 
sults for TV = 100. 

The energy of the next-to-leading order large TV and 
Hartree approximations is the same as the exact one, but 
the energy of the leading order large TV approximation 
differs from it by terms of order 1/N. (This is because 
we need to keep the initial values of the parameters the 
same.) To make a comparison between the approxima- 
tions this difference has to be compensated for; we do this 
by rescaling time by a constant multiplicative factor so as 
to match the last oscillation maxima. This effect is of or- 
der 1/N. For TV = 100, we find that for < t < 100, the 
next-to-leading order 1/N approximation is always more 
accurate than the leading order; however, when compar- 
ing the next-to-leading order with the Hartree, although 
less accurate for t < 50, the Hartree approximation starts 
becoming more accurate at t ~ 50 (however, the errors 
are very similar in magnitude, of the order of a few per 
cent). 

We now return to a discussion of the blow ups first 
encountered in the N — 20 case discussed above. The 
failure of truncation schemes (of which 1/N is an exam- 
ple) to maintain a rigid connection with the existence of 
a probability distribution function is a well-known prob- 
lem in nonequilibrium statistical mechanics. It is often 
the case that, when this connection is lost (failure of real- 
ity or positivity conditions), instability soon follows. In 
our case, the violation of unitarity is manifest in that 
positive expressions such as (r 2 )/TV can turn out to be 
negative. Note that since both the Hartree and leading 
order 1/N approximations are variational in nature, they 
can never violate unitarity. In Fig. |^, we show the blow 
up or failure time for the next-to-leading order large N 
approximation as a function of TV; the failure time being 
defined as the time at which (r 2 )/N becomes negative. 
It is interesting to note that for values of TV near the 
critical value TV C = 18.60 where the effective potential 
extends down to the origin, the failure time starts to in- 
crease rapidly (consistent with our interpretation of the 
connection of the static effective potential to an associ- 
ated quantum state). For values of N greater than 21, we 
could not find failure for t < 150. Thus, it may well be 
that above a certain value of N, the unitarity violation is 
pushed out to times of no practical significance, or may 
even disappear altogether. 
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FIG. 8. Failure time for the next-to-leading order large TV 
approximation, as a function of TV, for the same set of param- 
eters as in Fig. |^. Here failure time is defined as the time at 
which (r 2 )/N becomes negative. 

We discussed previously how to choose initial condi- 
tions so as to reach the large N limit in a controlled man- 
ner. Starting with these initial conditions, we show in 
Fig. ^ results from numerical solution of the Schrodinger 
equation for the time evolution of (r 2 )/N. In the strict 
large N limit one expects pure harmonic oscillations 
about the minimum of the infinite- N effective potential. 
For any finite value of N, however, at sufficiently large 
times the pure harmonic motion is overcome by nonlin- 
ear effects, and interesting behavior is found, as shown in 
Fig. [l^, indicating the presence of a very non-Gaussian 
wave function. The ability to capture this long time be- 
havior is an important test of 1/N methods. 
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FIG. 9. Exact solutions for (r 2 ) /N as a function of time. 
From top to bottom, N = 1, 5, 10, and 20. 
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FIG. 10. Very long-time behavior of the exact results for 
(r 2 ) /TV for < t < 500. The top figure is for N = 1 and the 
bottom figure for iV = 10. 

In order to test whether Hartree or the next-to-leading 
order approximation incorporate nonlinearities correctly 
so as to capture the late time modulation behavior, we 
ran a comparison against the numerical results for N — 
21, the results being displayed in Fig. (|TT|) . It is clear 
that both approximations do not give satisfactory results. 
This provides additional motivation for the development 
of alternative 1/N expansions which would incorporate 
selective resummations in order to reduce the coefficient 
of the error term at late times. 
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FIG. 11. (r 2 )/N for N = 21, for the exact, next-to-leading 
order large N, and Hartree approximations for late times. 
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B. Shifted Gaussian Initial Conditions 

We now discuss the time evolution of a qua ntum state 
having an initial wave function given by Eq. ( 7.16| ). For 
this problem, because of the lack of symmetry, exact so- 
lutions were only obtained for N < 2. For N = 1, de- 
pending on whether the energy is above or below the 
barrier height, one observes either slow tunneling with 
rapid oscillations in one well, or slower oscillation in the 
complete range. At these low values of N, the large-TV 
expansion breaks down quickly, as in the quantum roll, 
but even here at N = 1, the 1/N corrections improve the 
short time accuracy of q(t). 

A more relevant comparison is to consider larger values 
of N at which the approximations have a better chance 
of capturing the exact behavior. In the next figure, we 
compare the Hartree with the leading and next-to-leading 
order large- N approximation for q(t), at N — 50. Fig. [l2| 
displays the results for a run with E > £& using the 
equal time Green's function approximation (see Sec. 0) 
method for obtaining the Hartree results. Unlike the sit- 
uation for the roll initial condition where the Hartree and 
next-to-leading order large N results are not dramatically 
different, here the qualitative behavior is quite dissimilar 
(whereas the Hartree and leading order results are in fact 
very close). 
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FIG. 12. Plot of q(t) vs t for E > E b with N = 50. 



X. CONCLUSIONS 

Testing the 1 /N approximation in quantum mechanics 
has already enabled us to arrive at some useful conclu- 
sions. In order to interpret our results, it is important to 
keep in mind that 1/N approximations are a form of rc- 
summcd perturbation theory and are therefore only valid 
at weak coupling. Thus for couplings of order unity, it is 
unrealistic to expect the approximation to give good re- 
sults for small values of N. Our results have shown that 



at sufficiently large N the next-to-leading order approx- 
imation is a clear improvement over the leading order 
approximation, however, at late times this approxima- 
tion (as well as Hartree) fails to capture the nonlinear 
effects that lead to nontrivial amplitude modulation of 
the radial oscillations in the quantum roll problem. 

We have noted the presence of a finite-time breakdown 
in the evolution given by the next-to-leading order ap- 
proximation. This result is related to the fact that the 
large ./V expansion for the expectation values does not 
necessarily correspond to a positive semi-definite density 
matrix when truncated at any finite order in 1/N. (At 
lowest order, the 1/N approximation is equivalent to a 
Gaussian variational ansatz for the density matrix and 
does not have this problem.) This last aspect is already 
clear even in static situations such as the lack of a real ef- 
fective potential for all N in the next-to-leading order ap- 
proximation. This type of finite-time breakdown induced 
by unitarity/positivity violation has also been noted in 
simulations of quantum systems where the coupled equal 
times Green's function approach was truncated at fourth 
order p3] or where high order cumulant expansion meth- 
ods were used [ p3| . 

Two aspects of this breakdown deserve further men- 
tion: First, the time at which breakdown occurs appears 
to be strongly connected with the behavior of the effec- 
tive potential. For values of N not very much bigger than 
the critical value N c (beyond which the effective poten- 
tial exits over the entire range of y), the breakdown time 
increases extremely steeply and may even be pushed to 
times long enough to be no longer an obstacle to prac- 
tical calculations ( this still needs to be demonstrated). 
Second, it is important to point out that avoiding the 
breakdown via a partial resummation does not automat- 
ically guarantee better late time accuracy (or conver- 
gence) since such a scheme is also only next-to-leading 
order accurate. However, it will help in the sense that one 
may carry out simulations at smaller values of N, thus 
making it easier to compare against the late-time numer- 
ical solutions of the corresponding Schrodinger equation. 

One possible way of correcting the problem of a man- 
ifestly positive operator such as (r 2 ) becoming negative 
is to solve for the full Green's function Gij(t,t'): 

Gij{t,t') = G l3 {t,t') (10.1) 

-111 d< i / ^ G ik{t,h)^kl{tx,h)Qio{t2,t') , 
k,i Jc Jc 

rather than the next-to-leading order one as in (|3.7|). 
This equation is the exact equation one obtains by vary- 
ing the effective action and it contains terms of all orders 
in 1/N (thus, strictly speaking, one is no longer truncat- 
ing at some fixed order). 

However, just making this correction does not increase 
the time period during which the approximation is ac- 
curate. In order to extend the accuracy of the 1/N ap- 
proximation to late times, it appears necessary to use 
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a more robust approximation based on the Schwinger- 
Dyson equations. Several approximations of this sort are 
possible, which may both cure the positivity problem as 
well as lead to accurate results at late times. These will 
be discussed separately p4|. 
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